Nonlinear dual-core photonic crystal fiber couplers 
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We study nonlinear modes of dual-core photonic crystal fiber couplers made of a material with 
the focusing Kerr nonlinearity. We find numerically the profiles of symmetric, antisymmetric, and 
asymmetric nonlinear modes, and analyze all-optical switching based on instability of the symmetric 
mode. We also describe elliptic spatial solitons controlled by the waveguide boundaries. 
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Photonic crystal fibers (PCF) have attracted a lot of 
attention due to their intriguing properties, potential ap- 
plications, and the development of successful fabrication 
technologies [jj. PCFs are characterized by a conven- 
tional cylindric geometry with a two-dimensional lattice 
of air holes running parallel to the fiber axis. Such struc- 
tures share many properties of photonic crystals, associ- 
ated with the existence of the frequency gaps where the 
light transmission is suppressed due to Bragg scattering, 
as well as the guiding properties of conventional optical 
fibers, due to the presence of a core in the structure. 

Recent theoretical and experimental results reported 
the studies and fabrication of dual-core PCF structures 
for broadband directional coupling or polarization split- 
ting H H H H 0. In PCFs, light confinement is re- 
stricted to the core of the fiber and therefore nonlinear 
effects, such as light self-trapping and localization in the 
form of spatial optical solitons pj, become important. In 
particular, similar to two-dimensional nonlinear photonic 
crystals |8j, a PCF can support and stabilize both fun- 
damental and vortex spatial optical solitons |9|, [l0| . In 
a sharp contrast with an entirely homogeneous nonlinear 
Kerr medium where spatial solitons are unstable and may 
collapse, it was shown that the periodic structure of PCF 
can stabilize the otherwise unstable two-dimensional soli- 
tons. 

In this Letter, we make a further step forward in the 
study of nonlinear effects in the PCF geometry, and ana- 
lyze the existence and stability of nonlinear guided modes 
and spatial solitons in dual-core photonic crystal fiber 
couplers. The beam propagation and power-dependent 
switching in nonlinear directional couplers have been an- 
alyzed for the planar waveguide geometry Here we 
generalize those results for PCFs, as well as study the 
existence and stability of guided modes and elliptic spa- 
tial solitons controlled by the PCF holes. In particular, 
we find numerically the profiles of symmetric, antisym- 
metric, and asymmetric nonlinear modes, and analyze 
all-optical switching based on the mode instability. 

We consider a simple model of PCF that describes, at 
a given frequency, the spatial distribution of light in a 
nonlinear dielectric material with a triangular lattice of 
air holes with the radius r in a circular geometry. We as- 
sume that the PCF material possesses a nonlinear Kerr 
response, and there are two holes at the center filled by 
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FIG. 1: (a-c) Schematic of three designs of a nonlinear dual- 
core photonic crystal fiber. The coupler is created by two 
missing neighboring holes in a nonlinear material. 



the same material creating a nonlinear defect, as shown 
in Figs.^a-c). In the substrate material, the linear re- 
fractive index is n s , whereas inside the holes it is n a . 
In the nonlinear regime, the light distribution in PCF is 
described by the equation, 

BE 

-i-Q^ = + W(x, y)E + V(x, y)\E\ 2 E, (1) 

where W(x, y) = n a + (n s — n a )V(x, y), A± = d 2 /dx 2 + 
d 2 /dy 2 is the transverse Laplacian, E is the normalized 
electric field, and V(x, y) is an effective potential describ- 
ing two solid cores in the lattice of holes. We normalize 
V = 1 in the material, and V = in the holes. 

To find stationary nonlinear modes of PCF, we look for 
solutions in the form E(x,y,z) = u(x, y) exp(i/3z), and 
obtain the nonlinear eigenvalue problem, 



f3u = A±u + W(x,y)u + V(x, y)\u\ 2 u. 



(2) 



In order to find the solutions of Eq. (0 for nonlinear 
localized modes, we consider a rectangular domain of 
the (x, y) plane and apply a finite-difference scheme, tak- 
ing respectively N and M uniformly distributed samples, 
Xi, < i < N and yj, < j < M of the variables, as well 
as the corresponding samples for the stationary state, 
Uij = u(xi,Xj) and the potential Vij = V(xi,yj). Substi- 
tuting these variables into the model J2J), and imposing 
homogeneous boundary conditions in all four edges of 
the domain, we obtain an algebraic nonlinear problem of 
2 x N x M equations with the same number of unknowns 
, which is finally solved by means of a globally con- 
vergent Newton method. The presence of the external 
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linear potential given by two missing holes and the lat- 
tice of air-holes makes the system non-scalable and its 
radial symmetry broken. Another approach, that takes 
an advantage of the lattice periodicity, was developed re- 
cently by Ferrando et al. 

First, we consider a single missing hole and find nu- 
merically solutions for the PCF spatial solitons 9]. Im- 
portantly, these stationary solutions are not perfectly ra- 
dial, but they are stabilized by the PCF holes, in a sharp 
contrast with the unstable self-trapped beams in nonlin- 
ear focusing Kerr media. In order to demonstrate this 
feature, we follow the standard analysis of the soliton 
stability Q and analyze the soliton power as a function 
of the soliton propagation constant. A positive slope of 
this dependence indicates the soliton stability. 
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FIG. 2: (a-c) Light distribution between the two PCF cores 
in Fig. ^b) for three distinct nonlinear modes: symmetric 

(A) , antisymmetric (B), and asymmetric (C). (d) Transverse 
profiles of the nonlinear modes (f3 = 3.95). 

Next, we study the dual-core nonlinear PCFs shown 
in Figs.nja,b), and find the families of the spatially lo- 
calized modes-the so-called PCF spatial solitons-as a 
function of the mode propagation number /3. The corre- 
sponding solutions are similar for two cases of Figs^a,b), 
and they can be envisaged as the modes of the effective 
dual-core fiber generated by the combined effect of the 
dual-core PCF refractive index and the nonlinear index 
induced by the mode amplitude itself. We find that the 
nonlinear dual-core PCF supports three distinct nonlin- 
ear modes, symmetric mode (A), antisymmetric mode 

(B) , and asymmetric mode (C), as shown in Fig. [3 for 
the case of Fig. ^b). The corresponding spatial profiles 
of these modes are shown in Fig. |2(d) as cross-section 
cuts along the line y = 0. All the calculations are for 
n s = 5, n a = 0, r = 0.75 and A = 2. 

To demonstrate the power relation between the modes, 
we follow the standard analysis of the soliton stability [T^ 
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FIG. 3: Bifurcation diagram of the coupler modes. Insert 
shows an enlarged part near the bifurcation point. Asym- 
metric mode C bifurcates from symmetric mode A above a 
certain threshold in the mode power. 



and plot in Fig. [31 the soliton power as a function of the 
soliton propagation constant, for the three different fam- 
ilies. We notice that only two modes, symmetric and an- 
tisymmetric ones, may exist for low powers, whereas the 
asymmetric mode bifurcates from the symmetric mode 
at a ceratin threshold value of the mode power, above 
which the symmetric mode becomes unstable. 
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FIG. 4: Switching curve calculated for the PCF nonlinear 
coupler of Fig. ^b). Due to instability of the symmetric 
mode, the light launched into one core only does not switch 
to the second core but remains in the same core. 

In order to study the switching properties of the non- 
linear coupler, we carry out a series of numerical sim- 
ulations using the standard beam propagation method. 
First, we calculate the stationary mode of a single- 
core PCF with the same parameters as the coupler cor- 
responding to different powers. Then, this mode is 
launched into one of the cores of the dual-core coupler 
and, after some propagation, the output power is calcu- 
lated at the same core. The propagation distance is de- 
termined considering the same coupler operating in the 
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FIG. 5: Transformation of the elliptic guided modes to cir- 
cular solitons for the growing mode power. Three examples 
on the top mark three points of the eccentricity curve for 
f3 = 4.2, /3 = 6, and f3 — 15 respectively. Inset: the corre- 
sponding beam cross-sections at y — 0. 

linear regime and propagating the linear single-core mode 
launched to one of the coupler cores up to the point where 
all the energy is transferred to the second core. For the 
parameters used in our examples this point is reached at 
z = 21.23. The switching curve is obtained by varying 
the input power for the fixed propagation length, and it 
is plotted in Fig. 21 Similar to other types of nonlinear 
directional couplers, the input power is transferred com- 
pletely to the second core for low powers, while for it re- 
mains in the initial one for higher input powers. A change 
between both the regimes takes place in a relatively short 
range of powers which constitutes a threshold where the 
power switching is triggered. For very low powers the 
energy is completely transferred to the second core, al- 
though from the plot in Fig. 21 a residual amount seems 



to remain in the first one. This nonzero behavior of the 
curve close to the origin is explained by the overlapping 
of the field of the second core due to its proximity. 

Finally, we study the third case of the closely spaced 
holes, shown in Fig. dc). In this case, there exists no 
bifurcation to the asymmetric state and the fundamental 
mode itself is elliptic as shown in the inset of Fig. [SJ In 
the nonlinear case, this elliptic guided mode gives a birth 
to an elliptic spatial soliton controlled by the boundaries 
of the holes. In Fig. [51 we plot the mode eccentricity pa- 
rameter e = [1 — (wx/wy) 2 } 1 / 2 , where w x and w y are the 
mode axes (widths), and thus quantify the transforma- 
tion of the elliptic guided modes (A) to the elliptic soli- 
tons (B) and then to the radially symmetric solitons (C). 
These results resemble the transformation of the shape 
of nonlinear guided modes in planar waveguides [l3|. 

In conclusion, we have demonstrated that several types 
of two-dimensional spatial optical solitons can be sup- 
ported by a nonlinear dual-core PCF structure with the 
Kerr nonlinearity. We have analyzed numerically the ex- 
istence and stability of symmetric, antisymmetric, and 
asymmetric nonlinear modes demonstrating that periodic 
refractive index of PCF provides also an effective stabi- 
lization mechanism for these composite localized modes 
to exist in a nonlinear Kerr medium, in a sharp contrast 
with an entirely homogeneous nonlinear Kerr medium 
where spatial solitons are known to be unstable, and they 
undergo the collapse instability. We have studied all- 
optical switching in the nonlinear dual-core PCF coupler 
associated with the instability of the symmetric mode. 
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